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iShocks: X-ray binary jets with an internal shocks model 
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ABSTRACT 

In the following paper we present an internal shocks model, iShocks, for simulating 
a variety of relativistic jet scenarios; these scenarios can range from a single ejection 
event to an almost continuous jet, and are highly user configurable. Although the pri- 
mary focus in the following paper is black hole X-ray binary jets, the model is scale 
and source independent and could be used for supermassive black holes in active galac- 
tic nuclei or other flows such as jets from neutron stars. Discrete packets of plasma 
(or 'shells') are used to simulate the jet volume. A two-shell collision gives rise to an 
internal shock, which acts as an electron re-energization mechanism. Using a pseudo- 
random distribution of the shell properties, the results show how for the first time it is 
possible to reproduce a flat/inverted spectrum (associated with compact radio jets) in 
a conical jet whilst taking the adiabatic energy losses into account. Previous models 
have shown that electron re-acceleration is essential in order to obtain a flat spectrum 
from an adiabatic conical jet: multiple internal shocks prove to be efficient in providing 
this re-energization. We also show how the high frequency turnover/break in the spec- 
trum is correlated with the jet power, uj, oc , and the flat-spectrum synchrotron 



flux is correlated with the total jet power, F v oc L w 
agreement with previous analytical predictions. 



1.4 



Both the correlations are in 
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,£h I 1 INTRODUCTION 



X 



, In recent years there has been a lot of interest in probing 
the disc-jet connection in a variety of astrophysical objects. 
The mechanisms behind the jet formation are still not fully 
understood, leaving many open questions: the origin of a flat 
spectrum (a ~ when F v oc v a ) is one such question. The 
flat spectra have been ob served both in act ive galactic nuclei 
(AGN) (for a review see |Cawthorne!ll99ll ) and X-ray bina- 
ries (XRBs) (|Fenderll200ll ); in the XRBs it has been seen 
to ex tend from Radio to near infra-red l|Corbel fe Fender! 
2002). It is thought that this spectrum originates from the 
jet via t he partially self absorbed sy nchrotron emission. 

The iBlandford fe Konigll |l979j) model attempts to ex- 
plain how such a flat spectrum could arise. In their model, 
the jet is assumed to be conical with the magnetic field per- 
pendicular to the jet axis and frozen in plasma. For a given 
frequency, an increase in the jet volume causes a decrease in 
the plasma optical depth. The inner, denser, parts of jets are 
optically thick to lower frequencies (e.g. radio): the higher 
the energy density of the jet volume, the higher the op- 
tical depth (for a given frequency). The radio frequencies 
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therefore peak in the outer parts of the jet, while the infra- 
red p eak in the inner parts of t he jet. The one drawback 
of the IBlandford fe Konigll l| 19791 ) model is the artificial re- 
plenish ment of the adiab atic energy losses suffered by the jet 
plasma. iMarscherl (|l98(t) showed a model that takes electron 
energy los ses into account, but are unabl e to reproduce a flat 
spectrum; iHiellming fe Johnston! (|l988T ) presented a model 
where it is possible to obtain a flat spectrum under slowed 
expan sion for a conical jet. A comprehensive s tudy bv lKaisei] 
l|200rj| (more recently [Pe'er fe Casellal I2009T ) shows that if 
the adiabatic losses are not replenished then i t is impossible 
to obtain a flat spectrum from a conical jet. iKaiseij (120061 ) 
also show that using a non-conical jet volume, it is possi- 
ble to recover a flat spectrum even with energy losses: the 
jet geometry requires fine tuning to minimize the adiabatic 
losses normally associated with a conical jet, while allowing 
enough of a change in volume to drive the changes in optical 
depth. 

In the following paper, we present a model that ad- 
dresses the problem of electron re-energization via a large 
number of 'small' shell collisions. Although the main focus 
in this paper is on black hole XRB jets, the model is not re- 
stricted to such systems: the model is scale independent and 
AGN jet volumes can also be simulated. The first part of the 
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paper goes through the details of the model, outlining the 
physics and the techniques employed. The results section is 
split into different jet scenarios: single ejection, double ejec- 
tion (single collision/internal shock), and multiple ejections 
(multiple internal shocks). The results from these increas- 
ingly complex scenarios are used to demonstrate the model's 
capabilities, in addition to exploring the internal physics of 
the relativistic jets. 



2 THE MODEL 



Our model is based on the ISpada et all l|200ll ) internal 
shocks model for radio-loud quasar. Many modifications, 
however, have been carried out to make the model more flex- 
ible, and applicable to different scales and scenarios. In our 
model the jet is simulated using discrete packets of plasma 
or shells. For simplicity, only the jets at relatively large an- 
gle of sight are treated. Each shell represents the smallest 
emitting region and the resolution in the model is limited to 
the shell size. While the simulation is running, the jet can 
'grow' with the addition of shells at the base as the previ- 
ously added shells move further down the jet. If the time 
interval between consecutive shell injections is kept small, a 
continuous-jet approximation is achieved. The variations in 
shell injection time gap and velocity cause faster shells to 
catch up with slower ones, leading to collisions: the inter- 
nal shocks, discussed later, are a result of shell collisions. A 
schematic of the model set up is shown in figure [T] the two 
conical frusta shown represent the shells. 

2.1 Shell properties 

The shell volume is based on a conical frustum (cone opening 
angle = jet opening angle, ^). Asa shell moves down the jet, 
it can expand laterally as well as longitudinally (figure [1]). 
The adiabatic energy losses are a result of the work done by 
a shell in expanding; implicit assumptions are made about 
the pressure gradient across the jet boundary that would 
result in a conical jet. The emitting electron distribution is 
assumed to be power-law in nature; each shell contains its 
own distribution. The power-law distribution is of the form: 



N(E)dE = K,E~ p dE 



(1) 



where E — 7771c 2 is the electron energy, p is the power-law 
index and n is the normalization factor. If the total kinetic 
energy density of the electrons, Ek, is known then k can be 
calculated for the two cases of power-law index: p 7^ 2, and 
p = 2. When p / 2, we have (with the electron energy is 
expressed in terms of the Lorentz factor with mc 2 = 1): 



E k = 
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and for p = 2: 

E k = k {[ln("/ max ) - ln(7 min )] + hw - 7mL]} , 



(2) 



(3) 



where the subscripts max and min denote the upper and 
lower energy bounds for the electron distribution. The re- 
lations given in equations [2] and [3] can therefore be used 



to calculate the change in electron power-law distribution 
when there is a change in the total kinetic energy density, 
assuming the power-law index and 7 m i n are fixed, ymin value 
throughout the following work is set equal to unity, while the 
power-law index is assumed to be 2.1. The electron energy 
distribution upper limit, J m ax> is initially set to be 10 6 , but 
allowed to vary with the energy losses. 

A magnetic field is essential to give rise to the syn- 
chrotron radiation. In the shells, the magnetic field is as- 
sumed to be constantly tangled in the plasma, leading to 
an assumption that the magnetic field is isotropic; hence, 
treat ed like an ultra-relativistic gas l|Heinz fc Begelmanl 
2000). If the magnetic energy density (Eb) is given, the field 
(B) can be calculated: 



E B 



2u 



(4) 



where /jo is the magnetic permeability. 

Other shell properties include the bulk Lorentz factor, 
r, and the shell mass, M. If there is a variation in the F 
of different shells in the jet, then the faster inner shells are 
able to catch up with the slower outer ones, causing shell 
collisions; the shell collisions create internal shocks, which 
ultimately generate the internal energy. 



2.2 Internal shocks 

When two shells collide, a shock forms at the contact surface. 
Some of the steps involved in two-shell collision, and the 
subsequent merger, are shown in figure [5] The collision are 
considered to be inelastic. With many shells present inside 
the jet, first we need to calculate the next collision time 
between two shells: a collision is said to occur when the outer 
boundary of the inner shell, R° uter , comes in contact with 
the inner boundary of the outer shell, _R*™" er . The following 
relation can be used to calculate the time interval for two 
shell collision: 



dt co u 



R 



Tjouter 



O-i) 



(•>) 



where the subscripts j — 1, j denote two consecutive shells, 
/3 e is the shell longitudinal expansion velocity (along the jet 
axis), and f3 is the shell velocity (= \/F 2 — 1/r); dt is cal- 
culated for all the shells inside the jet and the minimum 
of these collision time gaps is selected i.e. the next earli- 
est collision. The shell longitudinal expansion (/3 e ) is due 
to any thermal energy the shell may have. We do not ex- 
plicitly model a thermal electron population; however, the 
expansion effects of having such a population are incorpo- 
rated in the model . The shell expansion velocity is given by 
(|Spadaet alj|200ll ): 



r 2 i~(pi3' 3 y 



(6) 



where /3 S = v' s /c corresponds to the sound velocity in the 
plasma (in shell co-moving frame), and: 



3 M 



(7) 



with E' th being the shell thermal energy and M is the shell 
mass. The prime denotes quantities in the shell rest frame. 
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Figure 1. An illustration of shells in our jet model. If the outer boundary of the inner shell, (j), contacts the inner boundary of the 
outer shell, (j-1), a collision is said to occur. The lateral expansion is due to jet opening angle; the longitudinal expansion is due to the 
shell walls expanding within the jet. The illustration is not to scale. 



The lPanaitescu fc M eszaros ( 1999) treatment of shock prop- 
agation is followed to work out the various quantities as- 
sociated with the shock itself. The shock propagation can 
be split into two shock-fronts. The forward shock travel- 
ling from the contact surface and through the outer shell 
(j — 1). The reverse shock travels through the inner shell 
(?) (injected after shell (j — 1)). Once the shock-front has 
passed through, the plasma/shell is considered to be shocked 
(S) and will have different physical properties compared to 
the unshocked plasma (see figure [2}. In one of the shell co- 
moving frames (shell rest-frame), the shock- front (SF) ve- 
locity can be calculated as: 



{v's - i)( 7 r' s + 1) 



(8) 



/3T' s [7(r' s -l) + l] ' 

7 is the adiabatic index and Y' s corresponds to the shocked 
plasma and is given by: 



v's = r m r(i - p m p) 



(9) 



where the subscript m denotes the merged shell proper- 
ties. The merged shell is formed once the shock-fronts have 
passed through both shells, leaving one combined shell. The 
merged shell mass is simply the linear combination of the 
two merging shells i.e. 



Mr, 



(10) 



while the merge d shell Lorentz factor is given by 
jSpadaet al.ll200ll ): 



r - 



VU) r U) + M(j-i) r (j-i) 



1/2 



(11) 



vM(j)/r(j) + Mu-ij/ro-i) , 

where /i — M + r//c 2 and r\ is the shell internal energy. As 
the shell collisions are considered inelastic, the generated 
internal energy is given by: 



E in = VU) + Vu-i) + P(j) c ( r (j) 

+ M(j-i) c2 ( r (j-i) - r m ) . 



(12) 



Once the quantities, outlined above, associated with 
the shocked plasma are calculated, we are able to deter- 
mine the new merged shell length. Due to the shock prop- 



agation through the plasma, the two shell lengths cannot 
be combined linearly to give the merged shell length. The 
shock propagation has a compression effect on the shell; the 
merged shell length is given by: 
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and the density, p, is: 

r m 7 r' s + 1 
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(13) 



(14) 



The Lorentz factor, T, corresponds to one of the two shells 
involved in the collision while F' s is given by equation [9] We 
do not consider the re-energization of the new merged shell 
to be instantaneous, but instead the energy is dissipated 
over a time period the shock-fronts would take to cross the 
inner and outer shells combined i.e. 



dtf 



dl 



+ 



dl 



Prs Pfs 



(15) 



where the subscripts ER, FS and RS denote energy release 
(time), forward shock and reverse shock respectively. Due 
to the limitations in the way the shocks are modelled, the 
merged shell is given the corresponding length (equationll3[) 
at the point of creation; only the energy release is delayed 
over time dtER- If the adiabatic energy losses are taken into 
account, then the merged shell will also be losing energy 
during the period dtER- 



2.3 Adiabatic losses 

For an expanding conical jet, the adiabatic energy losses 
need to be taken into account. These energy losses are 
due to the work done by the jet material while expanding 
(dll = — PdV, where U, P and V are energy, pressure and 
volume respectively). Using the relation, P = nkT, and as- 
suming that all the synchrotron emitting electrons are highly 
relativistic, we have: 



To 



(16) 
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(a) 







(b) 


RS 

j 






j-1 



(c) 



(d) 

Figure 2. An illustration of a two shell collision leading to a 
merger: (a) when the inner shell, j, comes in contact with the outer 
shell, j-1, (b) a shock starts to form; (c) the forward shock, FS, 
travels through the outer shell while the reverse shock, RS, travels 
through the inner shell; (d) once shock fronts have traversed the 
two shells, a new 'merged' shell is formed. 



with the adiabatic index, 7 = 3; the subscripts '0' denote 
quantities before the change in volume; the Lorentz factor 7 
corresponds to the power-law electrons accelerated to high 
energies due to the shock-front passing through the plasma. 
The adiabatic energy losses for the kinetic energy contained 
in the power-law electron distribution can therefore be cal- 
culated using equation [16] To calculate the change in the 
power-law normalization , k, for a change in volume, the 
following relation is used: 



V 

"="'" 1 Vo 



(17) 



Once the power-law normalization is, initially, calculated us- 
ing equation [2] or [3] the above relation can then be used to 
calculate the subsequent changes in the normalization. If 
the shell is involved in another collision then the distribu- 
tion is re-calculated completely. The change in the maxi- 
mum Lorentz factor of the electrons, 7 ma;c , as the shell vol- 
ume changes, can be calculated using equation [T5] The com- 
bined effect of varying j m ax and n is to effectively 'evolve' 
the power-law electron distribution. 

The changes in the magnetic energy density can be de- 
termined in a similar manner. If we assume that the mag- 
netic field is constantly entangled in the plasma and treat it 



as an ultra-relativistic gas, we can calculate the changes in 
magnetic pressure, P(B) using: 



PV 1 = PoVo 7 
therefore, 



P(B) = P (B) 
where the adiabatic index 7 = 4/3. 



(18) 



(19) 



2.4 Partially self absorbed synchrotron emission 

To model the sy nchrotron radia tion, we employ the treat- 
ment outlined in lLongairl (| 19941 ). With only the power-law 
electrons present, the synchrotron emission calculation is 
simplified; the synchrotron monochromatic intensity is given 
by: 



Jy 



(1- 



(20) 



The emission coe fficient, J v , and the absorption coefficient, 
Xv, are given bv lLongairl (j 19941 ). These coefficients are a 
function of the power-law normalization and the magnetic 
field, which in turn depends on the energy density of the 
shells; r is the shell radius (r„ = Xv r )\ the Doppler factor, 



8 T = 



[r(l T /3(cos0)]- 



(21) 



where, 6 is the jet viewing angle and '=p' corresponds to 
either an approaching component or a receding component 
of the jet. 

If the shell area is given by A and the distance to the jet 
is D, then the flux, &F V , from a single shell has the following 
form: 

^ = 4-^^(1-—)- (22) 

The above relations are used to self consistently calculate 
the synchrotron spectrum as it varies with the shell prop- 
erties. A more indepth treatment of the radiative transfer 
should take the relativistic effects, such as the relativistic 
aberration, into account as they may affect the overall spec- 
tral normalization. However, with simplicity in mind, these 
effects are not taken into account here: any errors due to 
this approximation can be minimized by only treating the 
jets at a large angle of sight. 

When a shell expands, its optical depth (with respect 
to a given frequency) changes; hence, a shell that is opti- 
cally thick to radio frequencies can become optically thin to 
them, as it moves down the jet. The optical depth of the 
neighbouring shells is not taken into account when calculat- 
ing the emission from a given shell; only the jets at relatively 
large angle sight from the viewer can be modelled. In other 
words, only the emission directly from each individual shell 
is modelled. 



2.5 Model parameters 

So far, only some of the physics and the principles behind 
the model have been outlined. It is therefore important to 
list some of the parameters used in our model to see how 
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they influence the physical properties of a jet (appendix |X] 
contains all the model parameters). 

The internal energy of the shell is split between the 
electron kinetic energy (u e ), shell thermal energy (u t h)and 
magnetic energy(us) i.e. 

E int — U e + U t h + UB ■ (23) 

Equipartition between the electron kinetic and magnetic en- 
ergy is assumed. The thermal energy density, u t h, is assumed 
to be solely responsible for the longitudinal expansion of the 
shells. Although the power-law electrons would also exert 
pressure for a similar effect, this becomes more important 
when u t h is zero. However, any realistic jet scenario is mod- 
elled with the thermal population present; the zero-thermal- 
energy scenario is used mainly for demonstrative purposes. 
Also if the shell mergers are not taking place and all the in- 
jected shells are identical, the longitudinal expansion needs 
to be suppressed, otherwise equation [TT] generates an erro- 
neous value for the merged F m . Of course, a more detailed 
treatment for any future work, take into account not only the 
power-law electron pressure, but also model the synchrotron 
emission from the thermal population. 

It maybe possible to constrain the distance to the source 
(D) and the jet viewing angle (9) from observations. The 
shell mass on the other hand offers a free parameter. In case 
of massive ejection event, the mass can be set manually. For 
a continuous jet (multiple ejections), we have a choice of set- 
ting the individual shells' properties manually or sampling 
from a pseudo-random distribution of these parameters. In 
the pseudo-random case, if the jet kinetic luminosity, Lw , is 
known then it can be used to generate the shell mass values 
by using: 



Y.MjT.c 2 =L w t 



(24) 

3=1 

where tj e t is the duration the jet is "on" . Therefore the to- 
tal relativistic mass of all the shells present in the jet has 
to correspond to the jet kinetic luminosity of the jet being 
modelled. The time gap,di n j, between any two shells can 
also be set either manually or sampled from a Gaussian dis- 
tribution with a given mean and standard deviation. The 
shell bulk Lorentz factor (r), if not set manually, is picked 
from a random distribution of values with the maximum 
and the minimum values (r woJ) r TO j n ) set by the user; the 
bulk Lorentz factor for a shell does not vary as it moves 
through the jet, unless it is involved in a collision. The shell 
length (dl) and the jet opening angle (ip) are not constrain- 
able from observations. However, at least in the case of the 
shell length, attempting to achieve a continuous jet approx- 
imation results in the following relation: 



dl — / scale dti n j f^sh&llC 



(25) 



where l sca ie is a scaling factor with a maximum of unity. 
hcaie = 1 means no spatial gap between two consecutively 
injected shells; in the simulation, however, hcaie < 1 in order 
to avoid a 'pile-up' of shells at the source. 



3 RESULTS 

With the aim of demonstrating a few of our model's capa- 
bilities, the following section is split into three main scenar- 
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Figure 3. Emission from a single shell: radio and infra-red fre- 
quencies are shown. The shell optical depth, t v , is also shown; 
the horizontal line marks optical depth of unity. 



ios: single ejection, double ejection, and multiple ejections. 
Within each scenario, the effects of adiabatic energy losses 
are also explored. As stated earlier, the focus in the following 
sections is on black hole XRB jets; accordingly various shell 
properties, outlined above, are calculated and estimated for 
XRB scaled jets. 



3.1 Single ejection 

A single shell ejection scenario can be used to simulate the 
massive ejection events observed in X-ray binaries, in the 
form of radio flares, when the source is transi tioning from a 
hard state to a soft state ([Fender et al.l[20o3 ). When shells 
are injected into the jet, they only start emitting radiation, 
or "light-up" , when they have been involved in a collision i.e. 
the shock front has passed through and energized the plasma 
within the shell (unless if they are injected with internal en- 
ergy). In the case of a single shell, however, no collisions can 
take place; thus, the shell will require lighting-up artificially. 
This is achieved by creating a "shock-zone" at an arbitrary 
point, x s hock, along the jet; the shock-zone picture is remi- 
niscent of the jet models involving a single-shock-zone in the 
jet (see lFalckelll99rl iFalcke fc Markofj 12001 iMarkoff et ail 
l200ll . lMarkoff et al-fool iPe'er fc Casellall2009l and the ref- 
erences therein). In our model, once a shell passes through 
the shock zone it is energized instantly. A fraction of the 
shell's relativistic energy, Ef rac , is used as a scaling for the 
amount of internal energy given to a shell after energization; 
the shell velocity, F, and mass, M, remain unchanged. 



3.1.1 Without energy losses 

Not taking the shell energy losses, as it expands, into account 
is similar to assuming that any energy l osses are continually 
replenished (jBlandford fc Konigll mm . With this kind of 
set up, a shell is allowed to propagate down the jet and 
expand laterally. The shell longitudinal expansion, due to 
the thermal energy of the plasma, is suppressed: the thermal 
energy is set to zero. The change in volume is therefore not 
associated with any work done by the shell. 

The emission at radio and infra-red frequencies, from a 
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Figure 4. Infra-red and radio lightcurves for a single shell with- 
out adiabatic energy losses. 



Table 1. The parameters used for single ejection scenarios. 



Parameter 


fig ei m 


fig ei m 


<P 


5° 


5° 





40° 


40° 


D 


2 kpc 


2 kpc 


M 


1 X 10 7 kg 


1 x 10 7 kg 


r 


2.0 


2.0 


di 


1 X 10 4 m 


1 X 10 4 m 


Me 


0.5 


0.5 


»C 


0.5 


0.5 


Uth 


0.0 


0.0 


shock 


0.2 Is 


0.2 Is 


frac 


0.3 


0.3 


Sim. Duration 


7 x 10 4 s 


500 s 



single shell, is shown in figure [3] The plot also shows how 
the shell optical depth for the two frequencies changes as the 
shell moves down the jet. The increase in shell volume causes 
the shell to become optically thinner to lower frequencies. 
This behaviour is evident in the figure, as the radio emission 
peak is much further down the jet than the infra-red peak: 
the emission peaks at r„ ~ 1. Various shell and jet volume 
parameters are outlined in table [T] 

The radio and the infra-red lightcurves, in figure [4] il- 
lustrate how the radio rise time is much longer than that of 
the infra-red. These rise times are determined by the shell 
energy density, and how quickly this energy density changes 
with the change in volume. Therefore, in the case shown 
here, the shell properties are such that the infra-red peaks 
rapidly after the injection, but the radio takes much longer. 
It should be noted that because energy losses are not consid- 
ered here, only the energy density of the shell, due to lateral 
expansion, is decreasing; the shell's total energy content re- 
mains fixed. 



3.1.2 With adiabatic energy losses 

When the adiabatic energy losses are taken into account, 
an increase in shell volume causes a decrease in the shell 
internal energy. In this particular simulation, as in the pre- 
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Figure 5. Emission from a single shell with adiabatic energy 
losses: radio and infra-red frequencies are shown. The shell optical 
depth, t„, is also shown; the horizontal line marks optical depth 
of unity. 
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Figure 6. Infra-red and radio lightcurves for a single shell with 
adiabatic energy losses. 



vious scenario, the shell thermal energy is again set to zero. 
Therefore any increase in the shell volume is purely due to 
the lateral expansion of the jet. 

The radio and infra-red spectra from a single shell, with 
the adiabatic energy losses taken into account, are shown in 
figure [5] When compared with the spectra with no energy 
losses, figure [3] there are two main differences that become 
apparent: the radio peak is at much lower flux value, and at 
much smaller distance along the jet. The infra-red peak flux 
values, however, remain relatively unchanged; the peak oc- 
curs at smaller radii when compared to the no energy losses 
case. These differences can be explained when one looks at 
how the shell energy density changes with and without the 
energy losses present. 

The shell optical depth is a function of the shell's energy 
density, i.e. a shell becomes optically thin to lower frequen- 
cies as the energy density drops. When energy losses are 
not considered, the change in the shell volume solely effects 
the change in the energy density. However, when the adia- 
batic energy losses are taken into account, the energy losses 
along with an increase in the shell volume drive the change 
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Table 2. The properties of the two shells injected in the double- 
ejection scenario shown in figures [JJ E] EJind \To\ 



Parameter 


shell 1 


shell 2 


Inj. time 


1 s 


10 s 


Mass 


1 x 10 6 kg 


1 X 10 10 kg 


r 


2.0 


3.0 


d/ 


1 X 10 5 m 


1 X 10 7 m 



in the shell energy density. In such a scenario, there is a 
two fold effect on the shell energy density: increase in the 
shell volume and the decrease in the shell internal energy. 
This then means that the shell optical depth, for a given 
frequency, changes more rapidly than when only the volume 
is effecting the change. The emission intensity, however, is 
also affected by the energy losses. With the adiabatic energy 
losses active a shell may be able to peak, for example, in ra- 
dio frequencies at smaller jet radii, x, but the peak intensity 
is lower due to the energy losses suffered by the emitting 
electrons and the magnetic field. (If longitudinal expansion 
is also taken into account then the energy losses are accel- 
erated and a shell is able to peak in radio frequencies even 
earlier with a further decrease in flux values). This is can 
be seen in figure [5] and the lightcurves in figure [6] where 
the peak radio flux is nearly two orders of magnitude lower 
when compared with the peak radio flux in figure [3] The 
infra-red peak flux on the other hand suffers a relatively a 
small reduction. This is due to the shell becoming optically 
thin to IR frequency very quickly in both cases, thus not 
having the time to suffer much energy energy losses. When 
the adiabatic losses are active, the shell starts off optically 
thin to IR. This is because the relative change in volume 
from the moment of injection to the subsequent time step 
being sufficient to drop the shell energy density below the 
limit for the shell becoming optically thin to infra-red fre- 
quencies. The initial volume of a shell, in addition to relative 
change in volume, also plays a part in how quickly that shell 
becomes optically thin to a given frequency: a large enough 
shell could start off as optically thin to infra-red frequencies. 
The sharp cut off for the infra-red flux, in figure 4, is due to 
"/max dropping below the energy threshold for the power-law 
electrons to emit in the infra-red. 




Figure 7. The infra-red (1 X 10 14 Hz) emission from the two 
injected shells that later on merge (at ~100s) to become a single 
shell. The properties of the two injected shells are outlined in 
table [2] while the simulation parameters are outlined in table [3] 
Adiabatic energy losses due to lateral expansion only are being 
modelled. 




Figure 8. The radio (1 X 10 9 Hz) emission from the two shells 
injected that later on merge (at ~100s) to become a single shell. 
The properties of the two shells are outlined in table [2] while the 
simulation parameters are outlined in table [3] Adiabatic energy 
losses due to lateral expansion only are being modelled. 



3.2 Double ejection 

A double ejection scenario involves two shells being injected 
into the jet volume with a time interval, dti n j , between them. 
If the second shell, j, has a higher velocity than the preced- 
ing shell, j — 1, then the two shells should eventually collide. 
This scenario also demonstrates the core principle of mul- 
tiple ejections: a large number of two-shell collisions taking 
place all along the jet to give rise to multiple shocks. 

The shock-zone location, mentioned in the previous sce- 
nario, is set at zero. This means that the shells are injected 
with some internal energy instead of gaining it after pass- 
ing through an arbitrary point. The adiabatic losses are also 
modelled, but only due to lateral expansion of the jet; the 
shell thermal energy is set to be zero. 

The properties of the two shells injected into the jet 
volume are outlined in table [2] The first shell, (j — 1), is not 



Table 3. The parameters used for double-ejection scenario. 



Parameter 
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Figure 9. The lightcurve for the radio emission. The two arrows 
on the left signify injection of the two shells; the third arrow shows 
the time of merger. Initially (after the second shell injection at 
10s) the lightcurve comprises emission from both the shells; later, 
after merger at ~ 100s, only a single shell exists in the jet. The 
properties of the two shells are outlined in table [2] the simulation 
parameters are outlined in table [3] The unusual iotji(time) is used 
for demonstrative purposes. 
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Figure 10. The lightcurve for the infra-red emission. The two 
arrows on the left signify injection of the two shells; the third 
arrow shows the time of merger. Initially (after the second shell 
injection at 10s) the lightcurve comprises emission from both the 
shells; later, after merger at ~ 100s, only a single shell exists in 
the jet. The properties of the two shells are outlined in table [2] 
the simulation parameters are outlined in table [3] The unusual 
log (time) is used for demonstrative purposes. 



only less massive than the following one, j, but also larger. 
The combination of these parameters means that (j — 1) 
becomes optically thin to lower frequencies sooner than (j). 
This can be seen in the figures [7] and [8] where the radio peak 
for shell ( j — 1) is at much smaller jet radii than that for shell 
j; in the case of IR emission, shell (j — 1) is already optically 
thin at those frequencies, when it is injected, while j reaches 
the peak IR flux later. The figures also show the point where 
the two-shell collision, or merger, takes place: it is marked 
by a sharp increase in both the radio and the IR emissions. 
The infra-red emission had in fact faded away completely 



by the time the merger took place; thus, demonstrating the 
re-energization aspect of these collisions/internal shocks. 

The lightcurves shown in figures [9] and \W\ illustrate 
the lag between the high and the low frequency peaks, al- 
ready seen from individual shells. In this case, however, the 
lightcurves show the total emission from the entire jet (only 
two shells for this scenario; multiple shells case is presented 
below). The radio emission not only lags behind the infra- 
red emission, but also has a much lower peak flux value 
(due to energy losses). Also, the radio emission rise and de- 
cay times are much longer than the infra-red times. In the 
case of infra-red, there is a sharp rise in the emission at the 
point of shell merger. As mentioned earlier, the infra-red 
fades away almost completely by the time the shell collision 
takes place. The merger, however, has a compression effect, 
thus increasing the energy density of the newly formed shell, 
causing it to start emitting in infra-red. The internal energy 
generated at the collision is still not sufficient to make the 
shell optically thick to infra-red; therefore, we do not see a 
slow rise in the infra-red flux from the merged shell. The 
picture is slightly different for the radio: the merged shell 
has high enough energy density that it becomes optically 
thick to the radio frequencies, leading to a slow rise in the 
flux. The other point to note, for the radio, is that after the 
second shell is injected into the jet, it takes a long time be- 
fore the maximum flux is reached. This is due to the second 
shell taking a long time to become optically thin to the radio 
frequencies. 

3.3 Multiple ejections 

The multiple ejections set-up is also demonstrated with and 
without adiabatic energy losses taking place. When the en- 
ergy losses are absent, the internal shocks/collision are also 
omitted as no energy replenishment is required; the colli- 
sions, however, are modelled when energy losses are incorpo- 
rated. The shell properties in the case of collisions set-up are 
sampled from a pseudo-random distribution (as described in 
section I2,5[l . 

3.3.1 Without energy losses 

If the time gap between ejection is small then an almost con- 
tinuous jet can be approximated by multiple ejections. All 
the shells are injected with the same properties (time gap, 
velocity etc.), hence no collisions take place. The simulation 
parameters are given in table [3] 

The spectrum shown in figure [TT1 shows how a flat spec- 
trum (for a specific frequency range) can be recovered when 
no energy losses are considered. This set up can be com- 
pared to a situation where one assumes a constant replen- 
ishment of any energy los ses by an unknown mechanism 
(Blandfor d fc Konigll Il979l ). The evolution of a radio and 
an IR frequency along the jet can be seen in figures [T2l and 
1131 respectively. The two frequencies show a very different 
behaviour: IR spectrum shows a constant decline while the 
radio spectrum peaks much further down the jet. A look at 
the optical depths for the two frequencies shows that the 
injected shells, and ultimately the jet, is optically thin to 
infra-red (j v << 1); in the case of the radio frequency the jet 
becomes optically thin at a large distance from the source. 
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Figure 11. Spectrum of a jet modelled using multiple ejections 
(without energy losses). An almost flat spectrum is achieved over 
a large frequency range. It should be noted that the low frequency 
turn over point is dictated by the duration of the simulation: 
longer simulation means a longer jet which means a lower fre- 
quency turnover, (see table [4] for parameters). 
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Figure 12. The infra-red emission of a multiple ejections jet 
(without energy losses). Optical depth corresponding to the infra- 
red frequency is also plotted (long dashed line). The wriggles at 
the end of the IR spectrum are a numerical artefact. (see table [4] 
for parameters). 



Table 4. The parameters used to demonstrate multiple ejections 
(without energy losses). 
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(The same behaviour was observed in the case of single ejec- 
tions as well) . The emission for a range of frequencies (radio 
< v < infra-red) are shown in figures [141 The different emis- 
sions shown in figure [141 follow the R T asi oc v~ x relat ion (as 
predicted analytically bv lBlandford Sz Konigll l|l979l )'). 

3.3.2 With adiabatic energy losses 

The internal shocks are a possible way to address the prob- 
lem of replenishing the energy losses in a jet. In the simula- 
tions presented below the shells are expanding both longi- 
tudinally and laterally; therefore the adiabatic losses can be 
extremely fast, making the flat spectrum difficult to obtain. 
The spectra from simulations where shells are not injected 
with any internal energy can be compared with the spectra 
from the simulation where shells are injected with internal 
energy: internal shocks are the only source of the internal 
energy production in the former case, whereas in the latter 
they serve to replenish the energy losses only. 
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Figure 13. The radio emission of a multiple ejections jet (without 
energy losses). Optical depth corresponding to the radio frequency 
is also plotted (long dashed line); short dashed line shows where 
the optical depth ~ l.(see table|4]for parameters). 

The spectrum shown in figure [TS] shows a highly in- 
verted spectrum, when the shells are not injected with any 
internal energy. The internal shocks taking place are not 
sufficient to produce the internal energy in addition to re- 
plenish the energy losses that are taking place, resulting in 
an inverted spectrum. 

On the other hand, the spectra given in figure [T5] illus- 
trate how it is possible to obtain an inverted/flat spectra 
even with adiabatic energy losses taking place. In order to 
achieve this, it is necessary to inject the shells with inter- 
nal energy. In other words, when the internal shocks are 
used to produce the internal energy, plus replenish the adi- 
abatic losses, a flat spectrum is not obtainable; however, 
when the internal shocks are used for energy replenishment 
only, the flat/inverted spectrum is recovered. It should be 
noted that the role played by the injected internal energy 
(Ef ra c > 0.0) is somewhat more complicated as it also influ- 
ences the longitudinal expansion rate of the shells, thus how 
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Figure 14. The emission as a function of the jet radius for a 
range of frequencies in a multiple ejections jet (without energy 
losses), (see table l4l for parameters). 
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Figure 15. Time averaged spectrum (~ Iks) from a multiple 
ejection jet with adiabatic energy losses. The shells involved un- 
dergo lateral and longitudinal expansion. The shells are not in- 
jected with any internal energy, (see table [5] for parameters). 



Table 5. The parameters used for multiple ejections jet with 
adiabatic energy losses. 
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Figure 16. Time averaged inverted/flat spectra (~ 2ks) from 
multiple ejection jets with the adiabatic energy losses. The three 
spectra correspond to different jet kinetic luminosities. The shells 
involved undergo lateral and longitudinal expansion. The shells 
are injected with internal energy, creating a much different spec- 
trum from seen in figure [TBI ( see table [5] for parameters). 
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Figure 17. Time averaged inverted spectrum (~ 2ks) from mul- 
tiple ejection jets with the adiabatic energy losses. The adiabatic 
energy losses are due to the lateral expansion of the shells. See 
table [5] for parameters. 

quickly a shell becomes optically thin to different frequen- 
cies. When a shell is not injected with any internal energy 
(Ef r ac = 0.0), the longitudinal expansion does not begin 
until after the first merger; this means the colliding shells 
are quite small, causing higher frequencies to dominate the 
spectrum. Figure [T7] shows a highly inverted spectrum for 
a jet where the shells are injected with some internal en- 
ergy, but no shocks/mergers are taking place; the adiabatic 
energy losses are also being modelled (this is an almost iden- 
tical scenario to the one shown in figure [TT1 except for the 
additional adiabatic energy losses here) . It is clear that both 
the internal energy and the shocks re-acceleration are nec- 
essary to achieve a flat/inverted spectrum whilst taking the 
adiabatic energy losses into account. 

We can also note in figure \W\ that the flux is corre- 
lated with the jet kinetic luminosity. This is because Ej rac 
is scaled according to the relativistic energy of the shell, 
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which is related to the mass and the bulk Lorentz factor of 
the shell; the mass is dependent on the kinetic luminosity 
(see equation I24|l , which ultimately means that the higher 
jet kinetic luminosity creates shells with higher internal en- 
ergy, thus producing greater flux. Higher energy density also 
means that a shell would be optically thick to higher fre- 
quencies. The effects of jet kinetic luminosity on the flux 
(and the spectrum) are degenerate with Ef rac parameter: 
a lower luminosity jet, but with the higher Ef rac value can 
produce similar results. This degeneracy extends to any pa- 
rameter that influences the internal energy of the shell at 
injection; for instance, the jet opening angle and the shell 
length (Iscaie) will also influence the form of spectra ob- 
tained. The spectra show in figure [16] conform approximately 
to the relation: F v oc L^ A . This is in agreem e nt wi th the 
relation found analytically bv lHeinz fc Sunvaevl ([2003) , stat- 
ing: F v oc I/w' 4 - 

The flat /inverted spectra produced have shown other 
interesting correlations: both the high and the low fre- 
quency turnover points in the spectrum correspond to cer- 
tain jet properties. In the case of the high frequency break, 
the higher the jet power (mainly Lw, but also Ef rac ), 
the higher the break frequency. The high frequency break 
scales approximately as: v\, oc L^' 6 , which is remark- 
ably close to previously observed and c alculated relation of 



v b oc Lw 7 jFalcke fc Biermannl [f995l : iMarkoff et all 120031 : 



iHeinz fc Sunvaevll2003l ) . The low frequency turnover, on the 
other hand, is also affected by the jet luminosity, but the 
re-energization by the internal shocks appears to play the 
biggest role: both the number of shells present in the jet and 
the collision radii of the shells influence the low frequency 
turnover. 



link the shell properties (such as time gap between ejections 
and the bulk Lorentz factor) with X-ray timing information 
for example. Using the X-ray lightcurves to drive such jets, 
we can then look at the infra-red lightcurves produced, we 
have an additional diagnostics for checking self consistency 
in the model. This investigation may also show how a "single 
blob" picture may arise, where the radio rise and decay times 
plus the flux) are very similar to the ones for the infra-red 
Mirabel et atll 19981 : iFender fc Poolevll2000l '). This is clearly 
not compatible with the single shell picture presented in this 
paper. It can then be assumed that it is not a simple scenario 
like a single large blob being responsible for the observed 
massive ejection events. We must therefore delve further to 
try and understand why this is currently not reproducible 
with our model. 

It is safe to conclude that the results presented in this 
paper do not exclude other re-acceleration models, but the 
internal shocks model holds much promise in being able to 
reproduce the often seen flat/inverted spectra in addition to 
opening up the avenue for studying jet timing properties. 
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4 CONCLUSIONS 

The results presented in this paper show how it is possible 
to reproduce a canonical flat spectrum even when a dis- 
cretized jet is used. We have also shown how flat/inverted 
spectrum is also reproducible if the internal shocks are used 
for energy replenishment. If the internal shocks are used for 
the initial electron acceleration, on top of replenishing the 
energy losses, the spectra become highly inverted (a > 0). 

The multiple internal shocks created by multiple ejec- 
tion into the jet volume can provide a considerable amount 
of energy to the shells. The results show that even with 
an essentially random distribution of shells velocities and 
injection times, adequate re-energization is possible: the 
flat/inverted spectrum is achieved. This is an important re- 
sult in furthering our understanding of the jet physics. We 
have also seen that the high frequency break in flat /inverted 
spectra is correlated with the jet power; the lower frequency 
turnover shows dependence on the number of shells present 
in the jet in addition to their collision radii. Further inves- 
tigation is required to quantify fully the relation that may 
exist between the jet properties and various break frequen- 
cies. However, the break frequency correlations seen thus 
far, are in agreement with the theoretical prediction as well 
as the observations. 

The results outlined above also hint at being able to tie 
the timing properties (X-ray to Infra-red in the case of X-ray 
binaries) with the jet physics. It should also be possible to 



APPENDIX A: THE CODE 

We are in principle open and receptive to others wishing (in 
collaboration) to utilize our code to model and test different 
jet scenarios. If you are interested, please contact the author. 

The following section outlines all the customizable 
parameters in our model. With efficiency and expandability 
in mind, the model is coded in C++. Effort has been made 
to minimize dependencies and use GNU software only. Once 
compiled, the code can read in all the parameters from a 
simulation parameters file and a shell parameters file; for 
any subsequent changes to the parameters, the code does 
not require re-compilation. 

The customizable parameters for the code are as fol- 
lows: 

Jet Luminosity: Used when shell properties are pseudo- 
random. This determines the shell mass based on how many 
shells need to be injected. [J/s] 

BLFjnax: Used when shell properties are pseudo- 
random. This sets the mean of a Gaussian distribution to 
be sampled from. [r max ] 

BLF_min: Used when shell properties are pseudo- 
random. This sets the mean of a Gaussian distribution to 
be sampled from. [F m i„] 
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shell_width_factor: Sets the initial shell size, along 
the jet axis using the relation outlined in equation 1251 [no 
units] 

jet_opening_angle: The full opening angle of the jet. 
[degrees] 

source_distance: Distance to the source being mod- 
elled, [kpc] 

EThermaLfrac: The fraction of the shell internal 
energy to be given to the thermal energy; causes the 
longitudinal expansion, [no units] 

EelecKin_frac: The fraction of the shell internal en- 
ergy to be given to the total electron kinetic energy; 
determines the power-law distribution parameters, [no 
units] 

EMagnet_frac: The fraction of the shell internal en- 
ergy to be given to the magnetic energy; affects the 
magnetic field strength, [no units] 

powerlaw jndex: The power-law index, p, of the electron 
power-law distribution, [no units] 

e_gamma_min: "/ m in of the power-law electrons. [7] 

e_gamma_max: "/max of the power-law electrons. [7] 

nu_min: v m i n for the frequencies being modelled; used 
when logarithmically spaced frequency range is used. [Hz] 

nujnax: v m ax for the frequencies being modelled; 
used when logarithmically spaced frequency range is used. 
[Hz] 

nu_points: Determines the number of points for the 
logarithmic frequency grid, [no units] 

individual_frequencies Switch to turn on logarith- 
mically spaced frequency range; takes min. and max from 
above, [take integer values: n=off (then uses two frequencies 
below; y=on] 

nu_l Used if only two frequencies being sampled. [Hz] 

nu_2 Used if only two frequencies being sampled. [Hz] 

increase_time_resolution A switch for increasing 
sampling, for radiative emission, of the jet at the time 
interval given below, [takes integer values: n=off (in this 
case the jet is sampled only at "events" ; y=on] 

step_resolution If the above switch is on, this deter- 
mines the sampling time interval, [s] 

total_run_duration Total run time of the simula- 
tion, [s] 

shell_inj_duration The length of time for the shells 
injection, [s] 



avg_ejection_gap Sets the mean of the Gaussian dis- 
tribution for sampling the time interval between shell 
injections, [s] 

use_shell_file A switch for reading a file with shell 
properties: injection time, shell mass, shell Lorentz factor, 
shell width, [takes integer values: n=off; y=on (when on, 
jet Juminosity, BLFjnax, shelLwidth, shelLinj .duration, 
and avg_ejection_gap are all deactivated)] 

shell_file Name of the file with shell parameters, [should 
contain 4 columns with the corresponding shell proper- 
tiesunjection time, shell mass, shell Lorentz factor, shell 
width.] 

write_results_file A switch to activate writing every 
time step to a file, [takes integer values: n=off; y=on] 

results_file Name of the file for the above switch. 

final_time_step A switch to activate writing final 
time step of the simulation, [takes integer values: n=off; 
y=on] 

lightcurve_file Name of the file for writing the light 
curve data. [Always written by default] 

in_vacuum_expansion A switch to deactivate adia- 
batic losses. 

inj_int_energy A switch to inject the shells with in- 
ternal energy, after they have passed the "shock location" . 

rel_mass_frac Scale the amount of internal energy given to 
the shell by a fraction of the of the shell's relativistic energy. 

shock_location The shock location. Used when in- 
jecting with internal energy, [light seconds] 

slow .energization A switch to activate slow ener- 
gization; the shells are not energized instantly, but given 
the energy over a length of time determined by the shock 
crossing time, [takes integer values: n=off, y=on] 
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